In [1]:
import pyart
from matplotlib import pyplot as plt
%matplotlib inline
In [4]:
filename = '/data/collection/sgpcsaprsurI7.00/sgpcsaprsurI7.00.20150506.091302.maint.cfrad.20150506_081905.871_CSAP_v1715_SUR.nc.maint'
radar = pyart.io.read(filename)
In [5]:
max_lat = 26.5
min_lat =25.4
min_lon = -81.
max_lon = -79.5
In [8]:
print(radar.fields.keys())
In [9]:
#create an instance of the class using our radar
display = pyart.graph.RadarMapDisplay(radar)
#create a Matplotlib figure
f = plt.figure(figsize = [17,8])
#now we are going to do a three panel plot, resolution is a basemap parameter and determines the resolution of
#the coastline.. here we set to intermediate or 'i' ('h' for high 'l' for low)
plt.subplot(2, 3, 1)
display.plot_ppi_map('differential_reflectivity',
vmin = -7, vmax = 7,
resolution = 'i')
plt.subplot(2, 3, 2)
display.plot_ppi_map('reflectivity',
vmin = -8, vmax = 64,
resolution = 'i')
plt.subplot(2, 3, 3)
display.plot_ppi_map('velocity',
vmin = -15, vmax = 15,
resolution = 'i')
plt.subplot(2, 3, 4)
display.plot_ppi_map('clutter_filtered_differential_reflectivity',
vmin = -7, vmax = 7,
resolution = 'i')
plt.subplot(2, 3, 5)
display.plot_ppi_map('clutter_filtered_reflectivity',
vmin = -8, vmax = 64,
resolution = 'i')
plt.subplot(2, 3, 6)
display.plot_ppi_map('clutter_filtered_copolar_correlation_coefficient',
vmin = 0.5, vmax = 1,
resolution = 'i')
In [ ]: